Abnormal number of Nambu-Goldstone bosons in the color-asymmetric 2SC phase of 

an NJL-type model 



D. BlaschkcE 

Fachbereich Physik, Universitdt Rostock, D-18051 Rostock, Germany 
Joint Institute for Nuclear Research, 141980 Dubna, Moscow Region, Russia 

D. EberlQ 

Institut fur Physik, Humboldt-Universitat zu Berlin, D-12489 Berlin, Germany 

K.G. Klimenkd 

Institute of High Energy Physics, 14.2281 Protvino, Moscow Region, Russia 
M. K. Volkov, and V. L. Yudiche\| 



O 
O 

(N 



Joint Institute for Nuclear Research, 141980 Dubna, Moscow Region, Russia 

We consider an extended Nambu-Jona-Lasinio model including both (qq)- and (gg)-interactions 
with two light-quark flavors in the presence of a single (quark density) chemical potential. In the 
color superconducting phase of the quark matter the color SU C {3) symmetry is spontaneously broken 
down to SUc(2). If the usual counting of Goldstone bosons would apply, five Nambu-Goldstone (NG) 
bosons corresponding to the five broken color generators should appear in the mass spectrum. Unlike 
that expectation, we find only three gapless diquark excitations of quark matter. One of them is an 
5'(7 c (2)-singlet, the remaining two form an 5 , (7 c (2)-(anti)doublet and have a quadratic dispersion law 
in the small momentum limit. These results are in agreement with the Nielsen-Chadha theorem, 
according to which NG-bosons in Lorentz-noninvariant systems, having a quadratic dispersion law, 
must be counted differently. The origin of the abnormal number of NG-bosons is shown to be 
related to a nonvanishing expectation value of the color charge operator Qg, reflecting the lack of 
color neutrality of the ground state. Finally, by requiring color neutrality, two massive diquarks 
are argued to become massless, resulting in a normal number of five NG-bosons with usual linear 
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dispersion laws. 

^ . PACS numbers: 11.30.Qc, 12.39.-x, 21.65.+f 

Keywords: Nambu-Jona-Lasinio model; Color superconductivity; Nambu-Goldstone bosons 

a: 



c3 



I. INTRODUCTION 



It is well known that, in accordance with the Goldstone theorem 0,0, TV Nambu-Goldstone (NG) bosons appear 
in Lorentz-invariant systems if an internal continuous symmetry group Q is spontaneously broken down to a subgroup 
H (here TV is the number of generators in the coset space Q/'H), i- e. the number of NG-modes is equal to the number 
of broken generators. However, in Lorentz-noninvariant systems the number of NG-bosons can be less than TV. In this 
case, the counting of NG-bosons is regulated by the Nielsen-Chadha (NC) theorem 0: Let n\ and be the numbers 
of gapless excitations that in the limit of long wavelengths have the dispersion laws E ~ |p| and E ~ \p\ 2 , respectively, 
then N < ni + Ini. (Here, E is the energy and p is the three-momentum of the particle.) In particular, this theorem 
is valid for relativistically covariant theories as well, since in this case: i) the total number of NG-bosons equals TV, 
the number of broken symmetry generators Q ; ii) evidently, the dispersion law for these TV massless excitations looks 
like E ~ \p\, thus TV = n x . 

Recently, in some relativistic models describing the dynamics of the kaon condensate in the color-flavor-locked phase 
of dense quark matter, an abnormal number of NG-bosons has been revealed 0,0- Since the Lorentz invariance is 
broken in this case and some of the gapless excitations have a quadratic dispersion law, there are no contradictions 
with neither Goldstone nor NC theorems. The superfluid 3 He in the A-phase and ferromagnets 0, H| are other 
known examples of condensed-matter systems with an abnormal number of NG-bosons. 

In the present paper, we demonstrate the abnormal number of NG-bosons in the dense color superconducting phase 
(2SC) of quark matter for a simple version of the Nambu-Jona-Lasinio (NJL) model with two light quarks and a 
single (quark number) chemical potential. In this phase, which can be realized naturally only at sufficiently large 
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values of the chemical potential (300 MeV< fi <1 GeV) , the initial color SU C (3) symmetry is spontaneously broken 
down to the SU C (2) group. Hence, in accordance with the usual counting of the Goldstone theorem, one might expect 
five NG-bosons, corresponding to the five broken symmetry generators, to appear in the NJL model. We shall prove 
that in the 2SC phase only three gapless excitations that can be identified with NG-bosons appear. Two of them 
have a quadratic dispersion law, thus yielding no contradiction with the NC theorem. For further applications, notice 
also the following important criterion which is sufficient for the equality between the number of NG-bosons and the 
number of broken generators If Q u i = 1, N is the full set of broken generators and if ([Qi, Qj}) = for any 
pair then the number of NG-bosons is equal to the number N of broken generators. 

Recall, that there is an alternative approach to study color superconductivity, which is based on the weak cou- 
pling perturbative QCD 9J. In this case, using the Schwinger-Dyson equation with one gluon exchange, the color 
superconducting phase was proved to exist at asymptotically high densities. In the two-flavored QCD investigations 
of color superconductivity five scalar NG-bosons are shown to exist, which are mixing with gluons and required as 
longitudinal components for five massive gluons (color Meissner effect) [lOj. Hence, the results concerning the number 
of NG-bosons of both the QCD and the NJL model approaches seem to be in contradiction 1 . The origin of the above 
discrepancy is related to the fact that the ground state of the considered simple NJL-type model is not color-neutral. 
Indeed, since the expectation value of the color charge operator Qs is non-vanishing, (Qs) ^ 0, the criterion [5j is 
not applicable, and an abnormal number of NG-bosons must arise. On the other hand, in quark models with a 
color-neutral ground state, according to the above criterion, there should arise a normal number of five NG-bosons 
all having normal (linear) dispersion laws. Most interestingly, we find from the analysis of the mass spectrum in 
our model that the masses of two scalar diquarks are proportional to the ground state expectation value of the color 
charge (Qs)- Thus, when color neutrality, i. e. (Q s ) = 0, is imposed in NJL-type models as an additional condition 
(which can be realized by introducing a color-chemical potential, fig, "by hand" [HI), then two additional massless 
particles should appear, and the total number of NG-bosons would b e eq ual to five. Recently, it was shown that the 
ground state of the 2SC phase of QCD is automatically color-neutral [lj due to a dynamical generation of the color 
chemical potential [is by gluon condensation. So, we see that the original discrepancy in the number of NG-bosons 
between the QCD approach and an extended NJL approach would disappear. 

The paper is organized as follows. In Section II we investigate the considered NJL-type model in the framework of 
the Nambu-Gor'kov formalism and derive an effective meson-diquark action for a finite chemical potential. Moreover, 
the gap equations for the quark and diquark condensates are numerically studied. Sections III and IV contain a detailed 
analysis of the massless and massive excitations in the diquark sectors and explicitly show the influence of the color 
properties of the ground state on the diquark mass spectrum. Section V contains a summary and discussions. Finally, 
dispersion laws for gapless diquarks are derived in Appendix. 

II. THE MODEL AND ITS EFFECTIVE ACTION 

It is well known that perturbative methods are not applicable in the low-energy/density QCD regions. Instead, 
effective field theories are usually considered. The most popular effective theories are based on Lagrangians with 
four-fermion interactions, like the NJL-type models. Let us first give several (very approximate) arguments somehow 
justifying the chosen structure of our QCD-motivated NJL model introduced below. For this aim, consider two-flavor 
QCD with a nonzero chemical potential and the color group SU C (3). By integrating in the generating functional of 
QCD over gluons and further "approximating" the nonperturbative gluon propagator by a S— function, one arrives 
at an effective local chiral four-quark interaction of the NJL type describing low-energy hadron physics. Finally, by 
performing a Fierz transformation of the interaction term and taking into account only scalar and pseudo-scalar (qq)- 
as well as scalar ((7g)-interaction channels, one obtains a four-fermionic model given by the following Lagrangian (in 
Minkowski space-time notation) 2 

L = q[j v id u + - m ]q + G^qqf + (qi^rq) 2 ] + G 2 [q C ; e s 1 5 r 2 q][qe s 1 5 r 2 q c ]. (I) 

In (Jl, /i > is the quark chemical potential, which in isotopically symmetric quark matter is the same for both 
quark flavors, q c = Cq 1 , q c = q l C are charge-conjugated spinors, and C = ij 2 j° is the charge conjugation matrix 
(the symbol t denotes the transposition operation). The quark field q = qi a is a flavor doublet and color triplet as 
well as a four-component Dirac spinor, where i — 1, 2; a — 1, 2, 3. (Latin and Greek indices refer to flavor and color 
indices, respectively; spinor indices are omitted.) Furthermore, we use the notations f = (r , r 2 , r 3 ) for Pauli matrices 
in flavor space; (e s ) a/3 = e a/3S is the totally antisymmetric tensor in color space, respectively. Clearly, Lagrangian 



Of course, it is necessary to point out that the QCD weak coupling considerations are valid at asymptotic values of the chemical 
potential, i. e. at fi > 1 GeV, whereas NJL model results are correct for sufficiently lower values fi < 1 GeV. Thus, a direct comparison 
of results seems to be problematic. 
2 The most general Fierz transformed four-fermion interaction includes additional vector and axial-vector (qq) as well as pseudo-scalar, 
vector and axial-vector-like (gq^-interactions. However, these terms are omitted here for simplicity. 
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||TJ is invariant under the chiral 5/7(2) rx SU(2)r (at m = 0) and color S U r {3) symmetry groups. The physics of 
light mesons [3 III El El , diquarks [H El and meson-baryon interactions |19|, |20| was successfully described in 
the framework of different NJL models. These effective theories were involved for the inve stig ation of both ordinary 
[21II22I l23l | . hot and dilute E!> an d color superconducting dense quark matter [23.l25ll2^.l27j| . Usually, on the basis 
of light-meson and baryon phenomenology, the following restrictions on the coupling constants are considered 

Gi > tt 2 /(6A 2 ), G 2 < Gx , (2) 

where A is the cutoff parameter in the three-dimensional momentum space, necessary to eliminate the ultraviolet 
divergences appearing when quantum effects (loops) are taken into account (usually A < 1 GeV). Since P is a 
low-energy effective theory for QCD, the external parameter /i is restricted from the above: /i < A. In the region of 
coupling constants J2Jl and at small values of \x (the case of low baryon densities), only the operator qq has a nonzero 
vacuum expectation value. Thus, the chiral symmetry of the model is spontaneously broken in this case, while SU C (3) 
remains intact. The behavior of quark and meson masses is established rather well in quark matter of low density 
|2lLl22| . In particular, the number of 7r-mesons, which are the three NG-bosons in this phase, is equal to the number of 
broken symmetry generators. Hence, this example is a good demonstration of the fact that Lorentz noninvariance is a 
necessary but not sufficient condition for the abnormal number of NG-bosons to appear in the system. At sufficiently 
high values of the chemical potential fi ~ 300-^350 MeV, the two-flavor color superconductivity phase occurs in model 
(1) [51. So, a nonzero diquark condensate ((qq) 0) is formed in this case, and, as a consequence, color symmetry 
is spontaneously broken down to the SU C (2) subgroup. A naive counting gives us five NG-bosons in this case (it is 
just the number of broken generators of the SU C (3) group). However, as it will be shown further, there are only three 
gapless bosonic excitations of the 2SC quark-matter ground state. Two of them satisfy the quadratic dispersion law, 
in agreement with the NC theorem. 

The linearized version of Lagrangian Q that contains auxiliary bosonic fields has the following form 



q\fid v + ^7° — a — m — 27'' 



Tir\q 



4Gi 



- + ^feWA] - ^VrV]. (3) 

Lagrangians Q and J2J arc equivalent on the equations of motion for bosonic fields, from which it follows that 

A s ~ q c e 5 j 5 q, a ~ qq, if ~ iqj 5 rq. (4) 

However, it is more convenient to start our consideration from Lagrangian J2J. Clearly, the a and if fields are 
color singlets. Besides, the (bosonic) diquark field A s is a color antitriplet and a (isoscalar) singlet under the chiral 
SU(2)l x SU(2)r group. Note further that a and A 5 are (Lorentz) scalars, but if are pseudoscalar fields. Hence, if 
(a) ^ 0, then the chiral symmetry of the model is spontaneously broken (at m = 0), whereas (A 5 ) ^ indicates the 
dynamical breaking of color symmetry. In the framework of the Nambu-Gor'kov formalism (for a recent relativistic 

treatment see, e.g. ref. |28|) quark fields are represented by a bispinor ^> = f g c)' Integrating in the generating 

functional based on the Lagrangian Q over the quark fields, one obtains an effective meson-diquark action of the 
original model (1) 



S cS (a,7f,A 5 ,A* 5 ) = - J d A x 



a 2 + if 2 A 5 A* 5 



4Gi 4G 2 



in , / D + K 



-Tr s/C2; m ( K + D - ) ■ (5) 



Besides of an evident trace over the two-dimensional Nambu-Gor'kov matrix, the Tr-operation in (JjjJ stands for calcu- 
lating the trace in spinor- (s), flavor- (/), color- (c) as well as four-dimensional coordinate- (x) spaces, correspondingly. 
We have used also the following notations: 

K+ = iA^eV , R- = -iA s e 5 j 5 , 

D ± =i'y v d v ±m Q -m -S ± , S ± = a ± i-fifr. (6) 

Let us introduce notations for the ground state expectation values of the meson and diquark fields: (a) = cr , (7?) = if Q , 
(A* s ) = A* 5 , (A s ) = A S Q and 



(Kt ) Dt ) Zt) = (K±,D±,i:±) 



<T=er ,7r=7f ,A <5 =A£,... 

The quantities er , ir*, A s a , A** 5 (k, S = 1, 2, 3) correspond to the global minimum of the effective potential V e g and can 
be found as a solution of the system of equations 

dVeS _ n dVcs _ „ dV c s _ dVcs _ „ 

dn k ' da ' dA s ' dA* s ' [ ' 
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usually called gap equations. In Q we used the following definition of the effective potential: 



o-,7f,A' 5 ,A' 5 =const 



-V cS (<j,tt,A s ,A s ) J d*x, (8) 



where, in the spirit of the mean-field approximation, all fields are considered to be independent of x. Without loosing 
any generality of the consideration, one can search for a solution of equations Q in the form: (A*, A 2 , A 3 )= (0, 0, A). 
Moreover, we suppose also that 7? = 0. (The last assumption is maintained by the observation that in the theory of 
strong interactions P-parity is a conserved quantity.) If A = 0, the color symmetry of the model remains intact. If 
A 7^ 0, then SU C (3) symmetry is spontaneously broken down to SU C (2), and the 2SC phase is realized in the model. 
In this case, the system of gap equations is reduced to the following one: 

(Jo f d 4 g f E+ E~ 2E+ IE' \ 

2G 1 ~ 1 J {2vyE\ q l-{E+y + ql-{E-y + D + {q ) + D_{ qo )}' (J) 
A f d^q f 1 11 



where D±(q ) = q 2 — (E ± ) 2 — | A| 2 , E ± — E ± fi, E = y q 2 + A/ 2 , and M — m + a a is the constituent quark mass. 

In these and other similar expressions, q$ is a shorthand notation for q + is ■ sgn(q'o), where the limit e — * + must 
be taken at the end of all calculations. This prescription correctly implements the role of /i as chemical potential and 
preserves the causality of the theory (see, e. g. 

Let us now make a field shift a(x) -> a(x) + a a , A* s (x) -> A* 5 (x) + A* 5 , A s {x) -> A s (x) + A s and then expand 

the expression © into a power series of the meson and diquark fields. The second-order term S^g of this expansion 
is responsible for the mass spectrum of mesons and diquarks. It looks like 



5^(a,7?,A 5 ,A* 5 )=- Jd A 



S A*S1 



A d A 



4Gi 4G 2 



where So is the quark propagator, which is an evident 2x2 matrix in the Nambu-Gor'kov space 

D+ K-\~ 1 fab 



whose matrix elements can be found by means of the projection operator technique [30l |: 

a= [jL3- e -W*-v){ ^- E+ 7 0A + + qo + E ~ 7°A 

J (2tt) 4 \ (? 2_ ( ^ +)2 + | A |2 e 3 e 3^ + ^ q 2 -(E-) 2 + \A\ 2 e 3 e^ 

b = -iAe 3 [ e - iq{x - y) { - <7 5 A + - 

J (2tt) 4 U 2 -(£+) 2 + |A|W 7 ~ + q 2 - (E-) 2 + |A| 2 e 3 , 

- - iA*e 3 f rf4g e -H(.*-v) J I , 1 

I ^-\4 C 1 „2 _ fJP+\2 , I A 12^3^3 ' ^ „2 _ ( j?-\2 



(2ir) 4 {q 2 ~(E+) 2 + \A\ 2 e 3 e 3 ' + q 2 - {E-) 2 + |A| 2 e 3 e 3 

d4 1 „-ig(*-y) f <lo + E+ oi . qo-E- - 

! ' ll - ~ _2 f-c-\2 , I A 12.3,3 ' J1 





(13) 




(14) 


3 7 5 A-}, 


(15) 


♦}■ 


(16) 



(2tt) 4 \q 2 - [E+) 2 + |A| 2 e 3 e 3 ' ~ q 2 - (E-) 2 + |A| 2 e 3 e 3 

where A± = | (l ± 7 B£zM2 ^ are projectors on the solutions of the Dirac equation with positive/negative energy, 

e 3 e 3 = diag(— 1, —1,0) is the projector (up to a sign) in the color space. Quantities (II 31) 11611 are nontrivial operators 
in the coordinate, spinor, and color spaces, but they are unit operators in flavor space. 

A tedious but straightforward analysis of shows that there arises a mixing between the a and A 3 fields. We 
have 

5^ ((t.tt, A s , A* 5 ) = 5$ (t?) A 3 , A* 3 ) +sf ) (A\ A* 1 ) + sf } (A 2 , A* 2 ), (17) 

where (p = 1,2) 

A* p ) = - / d 4 *^^- + Tp'fc* {aA"e" 7 5 rfA^eV} ■ (18) 
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In the following we are going to search for gapless bosonic excitations of the 2SC quark-matter ground state of the 
model only. Since we are mainly interested in the diquark sector, the expression for the 7r-meson effective action 

si 2 J (ir) is omitted here. Moreover, sf^ (a, A 3 , A* 3 ) has a rather cumbersome form, so we also do not show it. (If 
needed, both these expressions can be obtained immediately from Hill) .') 

Before proceeding further, we should fix the parameters of our model. Here, we choose the paramctrization procedure 
given in [lj] . We simplify that model by ignoring the corrections to the pion-quark coupling constant that come from 
the mixing between the pion and axial-vector. To fix the model parameters, we require the model to reproduce 
the Goldberger-Treiman relation Mg n — F n , where « 93 MeV is the pion weak-decay constant, and is the 
pion-quark coupling constant determined in the local NJL model as follows 



-iN r 



(2tt) 4 J {q 2 - M 2 ) 2 



(19) 



(This integral is divergent, and a regularization is supposed to be implemented here. In our work we use a three- 
dimensional cutoff to make such integrals meaningful.) 

Then, we require the observed pion mass (~ 140 MeV) to be reproduced by the local NJL model with our parameters. 
After this, one more condition is needed to complete the parameter fixing procedure. We consider two possibilities 
to finalize it: i) fix the model parameters so that the QCD sum rules estimate for the chiral condensate (qq) « 
(—245 MeV) 3 is reproduced in our model, ii) choose the model parameters in a way that allows one to obtain the 
correct description of the p — ► ittt decay, as it was done in 0] . In the first case (let us refer to this set of parameters as 
set- A), the cutoff A is about 618 MeV, the four-quark interaction constant G\ is equal to 5.86 GeV~ 2 . The remaining 
constant, G2, is chosen according to the Fierz-transformation as Gi — 3Gi/4 = 4.40 GeV ~ 2 . The current quark 
mass m is fixed via the gap equations to the value: 5.7 MeV. The constituent quark mass in the vacuum (at zeroth 
temperature and chemical potential) amounts to 350 MeV. 

In the second case, we obtain a different set of model parameters, which will be referred to, throughout the paper, 
as the set-B. Imposing the condition that an extension of this model to the vector channel would give the experimental 
value for the p-meson width (the p-meson decays mostly into a couple of charged pions, and the decay is described 
by the constant g p ~ 6.1.) and keeping in mind that in the local NJL model one gets a natural connection between 
g p and g^ (discarding pion-axial-vector transitions), V6g p — g^, we obtain A = 856 MeV, G\ — 2.49 GeV -2 , 
G2 = 3Gi/4 = 1.87 GeV~ 2 , and too = 3.6 MeV. In the vacuum the constituent quark mass is 233 MeV, which is 
small, compared to the case, where the chiral condensate is used to fix the model parameters 3 . A similar procedure 
has been implemented in V27\ in the chiral limit. Here, we keep the current quark mass non-zero in order to reproduce 
the pion mass in the vacuum. 

Solving the gap equations at a fixed chemical potential, we obtain the constituent quark mass M and the color gap 
A that satisfy the global minimum of the effective potential. In Figs. Q and [3 one can see the quark mass M and 
the color gap A vs. the chemical potential for parameter sets A and B. For the parameter set A, one observes a 
first-order phase transition from hadron matter to the 2SC phase. The gaps (constituent quark mass and color gap) 
change abruptly when the chemical potential reaches the vacuum value of M (see Fig. ^| . On the contrary, if the 
parameter set-B is chosen, the gaps reveal a rather smooth character near the phase transition, which is typical for a 
phase transition of the second-order (see Fig. EJ) . 



III. GAPLESS EXCITATIONS IN THE DIQUARK SECTORS 

Since A 1 in (|17|) is not mixed with other fields, let us, first of all, focus on the diquark A 1 -sector. If A 1 (a;) = 

(<P! (x) +itp 2 (x))/V2,A* 1 (x)^{fi (x) -iipz (x) ) /y/2, then the effective action s[ 2 ) ( A 1 , A* 1 ) from can be presented 
in the form: 

si 2) (A\A* 1 )=si 2) ( ( p 1 ,^) =~J d*xd*m(x)r«(a;-tf) W (tf). ( 2 0) 
where the 2x2 matrix T(x — y) is the inverse propagator of the (p\, (/92-fields (summation over k, I = 1, 2 is implied 

_____ (2) 

in 1201) ). Its matrix elements can be found via a second variation of S{ . 

r2o(2) 

T kl (x-y) = — 1 (21) 
d<pi(y)d<pk(x) 



A larger quark mass can be obtained if one takes into account the mixing between the pion and axial-vector. However, we omit it in 
our paper for simplicity and for the reason that this mixing is negligible in the 2SC phase. 
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FIG. 1: The constituent quark mass M (solid line) and the color gap A (dashed line) as functions of the chemical potential 
(parameter set- A). 




FIG. 2: The constituent quark mass M (solid line) and the color gap A (dashed line) as functions of the chemical potential 
(parameter set-B). 
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In momentum space, the Fourier-transformed components of the T-matrix have the following structure 

rirc?) = Ib(p) = 5 (IVa(p) + TvI(-p)) , r^(?) = -r^(p) = | (r^W - IVa(-p)) ■ (22) 

The derivation of these relations is given in Appendix, where the quantity Fa*a(p) is also presented (see eq. (|A5|) ) . 
The particle dispersion laws are defined by zeros of detr(p) in the po-plane, i- e. by the equation: 



detro) - r u (p)r 2 2(p) - r 12 { P )r 21 ( P ) = t a * a ( p )t a * a (- p ) = 0. (23) 

Clearly, the solutions of l(23() with positive/negative po correspond to particles/antiparticles. Let us first put P — 
(po, 0,0,0). In this case, detr(p ) is an even function of p , so the solutions of eq. (|2*3)) in the p^-plane might be 
searched for. They are the squared masses of excited states in the A 1 -sector of the model. Moreover, all these 
formulae are then greatly simplified. Indeed, it follows from (|A5|) at p 2 = that 



/d 11 II 
7^4 T-T T f+\ n I \ + 7 — T F -^ n , , \ = -4p ff(po), (24) 
(27T) 4 |_ (po + qo + E+)D + (qo) ( P o + qo-E )D-{q ) j 

where it is possible to integrate over qo, using the following prescription: (po + go) ^ (po + Qo) + ■ sgn(po + <7o) an d 
Qo ~~ * Qo + is ' sgn(go) with £ — * 0+ (see also comments after formula JTUjl ), and 

1 f d 3 q f 1 0(£T) e(-E~) 



H{Po) 2j (27rr\( PQ + E++E+)E+ + ( Po -E--E A )E A + ( Po -E-+E A )E A \- 

In l(25|) = y/ (E ± ) 2 + |A| 2 . Since the integral in the RHS of this equation is ultraviolet divergent, we regularize it, 
as the other divergent integrals, by using a three-dimensional cutoff A, i. e. q 2 < A 2 . Then, eq. ((23(1 is transformed 
to 

p 2 H( Po )H(-po) = 0, (26) 

from which it is evident that there is a solution Pq = 0, corresponding to a gapless excitation of the 2SC ground state. 
Since the function H(po) does not vanish at po = and /j, ^ 0, this solution may be identified in the case of finite p 2 
with a NG-boson with the quadratic dispersion law: po ~ p 2 /H(0) (see Appendix). Now, let us suppose that at some 
nonzero value po = —mi the function H(po) has a zero, i. e. H{— mi) = 0. Then, detr(po) has two zeroes po = ±mi, 
i. e. the point p 2 , = m\ is a solution of eq. 126(1 . and the second bosonic excitation of this sector has the nonzero mass 
mi. In the low-po expansion we have: H(pq) — H(0) + P qH'(0) + ■ ■ •, where 

1 f d 3 q / 6{E~) 6{-E~) 1 



m-\J 



(2»)» I (E- + (E- - El)Eg (B+ + £j )£+ 



2 J (2»)HBa(E- +Ei) s [E- - E^fE^ E+(£++E A . 



Thus, in this approximation mi = H(0)/H'(0). The mi vs. /i plots for particular values of the coupling constants 
G\, G 2 and the cutoff parameter A discussed at the end of section 2 (or, equivalently, for the values A and M from 
Figs-Handle are presented in Figs.|21andQ] In Fig.[3]the values of mi are calculated for /i = 350-^400 MeV when the 
parameter set- A is involved. The leftmost point corresponds to the condition of the phase transition (/i w 350 MeV). 
At smaller p the color gap equals zero and eqs. 1)27(1 cannot be applied, since they have been derived by using eq. 1(10(1 
with A ^ 0. Similarly, in Fig.0]the axis n begins with 250 MeV. At smaller fi, (but grater than 233 MeV), the value 
of mi is almost zero. Common to both parameter sets A and B is that the values of mi are small and thus justify 
the series expansion of the function H(p ) around po = 0. 4 Analogous results are obtained for the A 2 -sector of the 
model. There, as in the A 1 -sector, only one massless boson with the quadratic dispersion law as well as a massive 
one with the mass m 2 = mi are found. Evidently, the two massless and two massive bosons in the A 1 , A 2 -sectors 
form antidoublets with respect to the unbroken SU C (2) symmetry. 

Now let us consider the sector with mixing of the complex diquark A 3 -field and the scalar er-meson. Introducing 
new real fields: 4>(x) = (A 3 (ir)-|- A* 3 (x))/v / 2, ip(x) = i(A* 3 (x) — A 3 (x))/v / 2, we rewrite the effective action (here 
we omit the cumbersome expression for this quantity, however it can be easily reproduced from 1(111) ') in the following 
form: 

S%(a,A^A* 3 ) = -yd 4 xd 4 y(a(x),<t>(x)Mx)) IL(x-y) (a(y),4>(y)My))\ (28) 



4 The numerical investigation shows that the mass value rn\ found in the low-po approximation (sec Figs.|3]and^J differs from the exact 
solution of the equation H(mi) = by no more than f0%. 



30- 



„ 20 
> 



E 




10- 



-| , 1 , 

350 375 400 

\x [MeV] 



FIG. 3: The diquark mass mi = H(0)/H'(0) as a function of the chemical potential (parameter set-A). 
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FIG. 4: The diquark mass mi = H(0)/H'(0) as a function of the chemical potential (parameter set-B). 
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where H(x — y) is the inverse propagator of the cr-meson and the diquark fields <f>,i/). Clearly, it is a 3x3 matrix, 

whose elements can be obtained from (I28|) by taking all second variational derivatives of iS^g over the fields a, (j), and 
if>. We have calculated the momentum-space matrix H(p) (corresponding explicit formulae and relations are omitted), 
and then, putting p = (p , 0, 0, 0), found that det(II(po)) has only one zero at p\ = 0. Therefore, only one NG-boson, 
which is an 5i7 c (2)-singlet, exists in this sector (as in the A 1 , A 2 -sectors, each zero of det(II(po)) in the pp-plane 
means the mass squared of the ground-state excitation, this time in the mixed cr-A 3 sector). 5 

As it was mentioned above, in the 2SC phase of the model, where A ^ 0, the color SU C (3) symmetry of the ground 
state is spontaneously broken down to the SU C (2) group. Hence, in accordance with the Goldstone theorem, it is 
generally expected that five NG-bosons would appear in the theory. However, we have proved that in the A 1 , A 2 , A 3 - 
sectors of the model only three massless bosons exist, which can be identified with NG-bosons. Therefore, there seems 
to be a deficiency of two NG-bosons in the framework of the model under consideration. In spite of this fact, there are, 
however, no contradictions with the NC theorem, since two of the three found NG-bosons have quadratic dispersion 
laws for \p\ — ► 0. In addition, the model contains also light massive diquarks. 

Usually, in the framework of the NJL model, the 2SC phase is studied in the region of coupling constants (J2J) . In 
this case, at /z = we have a chirally noninvariant phase even if m — > 0, and the transition to the 2SC phase occurs 
at some finite value of the chemical potential. Formally, however, one can consider the region 

w = {(Gi,G 2 ) : G 2 > tt 2 /(4A 2 ), tt 2 (G 2 - Gi) > 2GiG 2 A 2 }, (29) 

where (even at fi = 0) the 2SC phase is realized [3J. As can easily be seen from (|27|) . in the case of /i = we have 
H(0) = and mi = rn 2 = 0, and the total number of gapless excitations (NG-bosons) in the A 1 - and A 2 -sectors of 
the model equals four. They form two massless 5'C/ c (2)-antidoublets. Apart from this, there is one NG-boson, which 
is an 5f/ c (2)-singlet, in the mixed cr-A 3 sector. Naturally, all these excitations have a linear dispersion law, as it 
should be in the relativistically invariant case (see Appendix). However, at an arbitrary small fi, particles from one of 
these antidoublets acquire nonzero masses, and the dispersion laws for the two massless particles from the remaining 
antidoublet are changed essentially from the linear to quadratic laws. 

Note finally, that in the above consideration we used the real/imaginary part parametrization for diquark fields: 

A 1 (x) = (<p 1 (x) + i<p 2 (x))/V2, A 2 {x) = (<fi(x)+i<p 2 (x))/V2, A 3 (x) - A = (</>(.x) + iip(x))/V2, (30) 

where A p (x) are the unshifted diquark fields from Lagrangian JJJ and A = (A 3 ). However, one can use also an 
alternative "polar coordinate" parametrization: 

( A 2 (z) Uexpj-zV^^)^ ( |=4=| ( 31 ) 

where A Q are Gell-Mann matrices and the summation over a = 4, .., 8 is implied (the form of an exponential multiplier 
in 131|) corresponds to the fact that A p (x) is an S'C/ c (3)-antitriplct). Comparing (|30|) and (|31|l at small £i,m we see 
that 

2 

<Pi = &) ¥2 = -£4, <Pi=&, <P2 = = ?7, ip = 7/|£8- (32) 

Further, one should insert H31[) into iJSJ and expand the resulting expression into a series of meson- and £ a ,rj fields. It 
is easily seen that in the second order of the new variables the effective action is a sum similar to l|17jl. In particular, 
it means that the diquark fields £4,^5 are decoupled from other fields, and their effective action looks like l|2(J|l with 
the replacement (|32f) used. So, some elements of the inverse propagator matrix for £4,^5 fields differ by a sign from 
^m{x — y) 12 ip. but the determinant of the whole momentum space matrix is not changed. Hence, the mass spectrum, 
the dispersion laws etc., in the £4, ^-fields are the same, as in the case with old fields (f\,(f2- Similar conclusions 
are valid for the rest of new variables, i. e. such physical characteristics of the model, as particle masses, dispersion 
laws etc., do not depend on the choice of field parametrizations, H30JI or H31|l. In particular, the number of abnormal 
NG-bosons, which equals three in the model under consideration, is a parametrization invariant quantity 6 . 



5 A detailed investigation of the meson/diquark spectrum including the cr-A 3 mixing is in preparation. 

6 In papers [HIH an abnormal number of NG-bosons in the same SU(2) X U(l) toy model was discovered. But in Q the parametrization 
I3UI was used for the scalar doublet, whereas in 5] a "polar coordinate" parametrization, similar to 1311 . was used. Evidently, their 
results coincide, i.e. the abnormal number of NG-bosons, particle masses, formulae for dispersion laws are the same. 
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IV. COLOR-ASYMMETRY AND MASSES OF SCALAR DIQUARKS 

As it was noted in [l^, the color superconducting phase is automatically color-neutral in QCD, but this does not 
hold for the NJL model under consideration. Massive scalar diquarks and, as a consequence, the abnormal number 
of NG-bosons in the model discussed in previous sections testify for this. As mentioned in the Introduction, for 
the number of NG-bosons to equal the number of broken generators of the underlying symmetry group (in our case 
SU C (3)), the ground-state expectation value for the commutator of any two broken symmetry generators Q a and Qi, 
must be zero: ([Q a , Qb\) — 0|- This criterion is not fulfilled for the NJL model considered here, namely, the ground 
state expectation value of the color charge operator Q$ — q r y°X 8 q does not vanish, and, therefore, the usual counting 
rule for NG-bosons is not applicable. 

Let us look at the function H(po) at po = 0. One can rewrite the expression (|27J) in a different form: 

, , 1 f d 3 q 1 / N E+ E~\ 

H{0) -2JWf^{- 2 ^- E) + ^-Wj- (33) 

This formula can be expressed in terms of quark densities. For this purpose, we need the expression of the thermo- 
dynamical potential of the system, which in the mean-field approximation is equal to the value of V e s in its global 
minimum point, supplied by the gap equations, i. e. 

n(p) = V cS (a ,7r ,A s ,A* s ), (34) 

where V e g was given in ©. Further, one can divide the thermodynamical potential into several terms: 

m = ±a M + Oi^ + ^ (35, 

c— 1 

where fi c (/-0 * s the contribution from the color-c quark, defined as follows: 

d 3 q 



W = -2 {\E+\ + \E-\). (36) 

The contribution from color-1 and color-2 quarks, ^i(2)(/-0 are equal due to the remaining SU C (2) color symmetry 
after the global SU C (3) symmetry group is spontaneously broken in the color superconducting phase. The third color 
quark, which does not condensate, gives fl 3 (p). Then, one can calculate the densities n c of color-1, -2, and -3 quarks 
separately without introducing additional chemical potentials for each color simply by differentiating £l c (p) over //, 
i. e. 



9/i 



1,2,3; (37) 



"2 -2 I -e-l^--^). "■: =4 / -^9(fx-E). (38) 



f d 3 q 






1 (2^) 3 







Eq. I|33|) can thus be transformed to 



(2tt) 3 



#(°) = (»i + «a - 2n 3 ) = ^g^l . (39) 



From eq. Q39J) and the fact that the diquark mass mi = H(G)/H'(0) is found to be nonzero (see Figs. 3 and 4), one 
concludes that the 2SC-ground state in the NJL model under consideration has finite color charge (Qs). Evidently, 
this fact is related to the inequality of color quark densities > for paired (1,2) and unpaired (3) quarks, i. e. 
to color-asymmetry. Moreover, as discussed in Appendix, the nonzero value of H (0) is responsible for the abnormal 
dispersion law p ~ p 2 /H(0) for two NG-bosons. 

To restore the local color neutrality in the NJL model, one can, e. g., include an additional term, fJ.sQ&, and impose 
the color neutrality condition: 

<Q 8 } = -^^0, (40) 
dps 

or consider the superconducting matter as composed of colored domains, with the global color charge being equal 
zero. The latter, however, goes beyond the mean-field approximation exploited here and should be considered as an 
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approximation for the case of large colored domains. Anyway, the true ground state should be checked to give the 
absolute minimum of the thermodynamical potential. 

The inclusion of an additional chemical potential "by hand", however, does not contradict QCD. As discussed 
in |12fl . a non- vanishing expectation value of the gluon field Aq, corresponding to the eighth generator in the SU C (3) 
group, gives rise to a non- vanishing value of /is, which then cancels certain tadpole contributions and restores the 
color neutrality. This cannot occur in the NJL model automatically because there are no gluons in it. 

V. SUMMARY AND DISCUSSIONS 

Recently, it has been shown that in QCD with a nonzero strangeness chemical potential the number of NG-bosons 
can be less than the number of broken symmetry generators 0,TE|, so that the usual counting rule for NG-bosons 
evidently does not hold. In the present paper, we have presented another relativistic model providing us with an 
abnormal number of NG-bosons. This is an NJL-type model with two light-quark flavors, where a single quark 
chemical potential is taken into account. In the color superconducting phase of this model, the color SU C (3) group is 
spontaneously broken down to SU C (2), i. e. the number of broken symmetry generators equals five. Despite of this 
fact, there appeared only three NG-bosons in the 2SC phase; two of them form an antidoublet with respect to the 
unbroken SU C (2) subgroup and have quadratic dispersion laws. The remaining one is an 5f7 c (2)-singlet, so there are 
no contradictions with the NC theorem [3j. Moreover, there exists an 5 l J7 c (2)-antidoublet of light diquarks A 1 , A 2 
with masses mi, the behavior of which at final chemical potential \i is shown in Figs. 3 and 4 for the parameter sets 
A and B, respectively. 

Most interestingly, we found that the diquark masses mi = H(0) / H'(0) are proportional to the ground state average 
of the color charge operator (Q$) = l/-\/3(ni +n 2 — 2n 3 ). Here, the quark densities n\, ni of paired quarks with color- 
1,2 are equal due to the remaining SU C (2) symmetry and larger than the density of the unpaired color-3 quark. 
Thus, the ground state of the considered NJL-type model is evidently color-asymmetric. It is just the appearance of 
this nonvanishing expectation value (Q%) which simultaneously leads to a violation of the criterion for the equality 
of the number of NG-bosons and broken symmetry generators Q , and to the appearance of the quadratic dispersion 
law for the NG-diquarks (compare eqs. (38), (A6)). 

In the investigations of color superconductivity in the two-flavored QCD, five scalar NG-bosons are argued to exist. 
They are mixing with gluons and required as longitudinal components for five massive gluons (color Meissner effect) 
[To| . The fact that we find an abnormal number of three NG-bosons does not directly contradict QCD. Namely, in the 
color superconducting phase of QCD a nonvanishing condensate of the eighth gluon field component A® appears, which 
induces a color chemical potential fis and cancels certain tadpole contributions, responsible for the nonvanishing color 
charge of the ground state [l^- This mechanism just leads to color neutrality. Obviously, there arises the question, 
how can one reconcile the considered NJL approach with QCD? Insofar as NJL-type models do not contain gluon 
fields, the required contribution from the condensed eighth gluon field does not follow automatically. Therefore, the 
tadpole contribution becomes unbalanced by gluon contributions. It can be be cancelled "by hand", e. g. with the 
help of the additional color chemical potential fig, simulating the omitted gluon contribution. In such enlarged NJL 
models the condition of color neutrality of the ground state, i. e. (Qg) = 0, can be imposed as an additional physical 
requirement ^lj. According to our results, one gets in this case two additional NG-bosons as well as a change of 
quadratic dispersion laws into normal linear ones. Thus, the abnormal number of three NG-bosons might, in principle, 
become converted into the normal number five found in QCD. 

Nevertheless, what to do now with these five NG-bosons arising in the color neutral 2SC-phase of an enlarged 
NJL-type model? Clearly, in the framework of the standard NJL model one cannot apply the color Meissner effect, 
in order to absorb these NG-bosons into longitudinal degrees of freedom of massive gluons. A possible way out of 
this problem could be an extension of the NJL mode via the inclusion of perturbative gluons, fluctuating around the 
original low-energy gluonic fields (realized e. g. by instantons or other non-perturbative background fields, which were 
already integrated out to yield the effective four-quark interactions) j^. Obviously, the color Meissner effect could 
then be realized for such perturbative gluons. 

In addition, let us remark that the above discussions refer to the case of locally neutral quark matter, where 
color neutrality should indeed be required. On the other hand, global color neutrality does not demand local color 
neutrality: for example, one can assume the ground state as being composed of colored domains. Anyway, to decide 
which of the considered ground states is in favor, the absolute minimum of the thermodynamic potential should be 
calculated in each case. (Note that the presence of domains means a violation of translational invariance and is not 
consistent with the simple mean-field approach.) 

Finally, notice that in the framework of the considered NJL-type of model there arises a mixing between the cr- 
meson and A 3 -diquark in the 2SC phase. The behavior of a-, A 3 - as well as 7r-meson masses vs. fi in the 2SC phase 
is now under consideration. 
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APPENDIX A: THE T-MATRIX ELEMENTS AND DISPERSION LAWS 



Using the effective action Ijl8(l . let us introduce auxiliary quantities 



Ta«a(z) 



2c(2) 



5 2 S. 



AA 



M 

SA 1 {y)SA* 1 (x) 4G 2 

S 2 S{ 2) _ 6{z) 
SA* 1 (y)5A 1 (x) ~ 4G 2 



iTr sc {a(z)e 1 ^d(-z)e 1 ^ 5 } 
iTT sc {d(z)e 1 1 5 a(~z)e 1 1 5 } 



(Al) 



where S(z) is the Dirac ^-function, z — x — y and a(z),d(z) are operators 113|) and l|16|) . respectively. Since 
8 2 s[ 2) /(SA 1 5A 1 ) = 5 2 s[ 2) /(SA* 1 8A* 1 ) = and 

5 S 



S<Pi V2\SA 6 A* J ' S(p 2 y/2 \5A SA* 
the matrix elements of the inverse propagator matrix 1|21J) have the form: 

r n {z) = r 22 ( z ) = ^(r A » A (2) + r AA .(z)), r 12 (z) = -r 21 (z) ^(r 4 . A ( z )-r AA .(z)). 

For an arbitrary function F(z), it is possible to define the Fourier-transformed one, F(jp), by the relation 

d 4 p 



F(p) 



(TzF(z)e v - 



F(z) 



(27T)' 



■F(p)e- 



(A2) 



(A3) 



Now, using l|A3|l . it is possible to get from HAlfl 



4G 2 



iTr, 



r A » A (p) 



q {a(q+ P )e 1 1 5 d{q)e 1 ^}, 
{d(q+p)e 1 1 5 a(q)e 1 1 5 }, 



(27T) 4 

d 4 q { - 

WT 4 



(A4) 



where the Fourier-transformed expressions d(q), d(q) can be easily derived from (|13fl . (|16|) . It follows from (|A4() that 
Fa*a(~p) = r A A*(p)- So, taking into account (|A2|) . one can easily obtain the relations l|22l) . 
After tedious calculations we arrive at the expression 



Ta*a(p) 



4? 



d 4 q 



4»Po 

(2tt) 4 7 [(po + ^ + zi) 2 -^] 



Po + qa+ H + E q p + q + f-t- E q 



D-(q ) 



D+(q ) 



2*) 



[(Po + 90 + p) 2 



s 9 p 2 + (go + /i + s 9 )(p-(j) s 9 p 2 -( go + /i-s 9 )(p-g) 



E q D_(q ) 



E q D + {q ) 



(A5) 



Here, the notations that were introduced after (jl()|l are used. Moreover, E q = E = \/ q 2 + M 2 . Note, that expression 
(|A5I) is obtained under the assumption that A ^ 0, i. e. it is valid only in the 2SC phase. In this case, the gap 
equation (|10f) was used to eliminate the coupling constant G 2 from relation (|A4I) in favour of other model parameters. 

Quasiparticle dispersion laws are defined by (|23|l or, equivalently, by the two equations r A » A (±p) = 0. One can 
easily see that at p 2 — the first integral in RHS of (|A5(I coincides with 124(1 having an evident zero at po = 0. This 
is just the dispersion law of a gapless particle (NG-boson) at p 2 = 0. Let us find the dispersion law for the NG-boson 
at small nonzero values of p 2 , i. e. try to solve the equation T A * A (p) = at p — > 0. Expanding expression (|A5|) 



into a Taylor-series of (po,P) at the point (0,0) and taking into account only leading terms, we have, at ^ ^ 0, the 
following equation connecting the energy and momentum of the massless particle (in section III, it was proved that 
such a particle, NG-boson, exists in the A 1 -sector of the model): 



r A * A (p) = 4ip H(0) 
Aip 2 



Aip 



,72 



d 4 q 



(2*)* J D^q )[(q + ^ 2 -E 2 } 
d 4 q f 2q 2 



2<f 



(2tt)4 J D + (q Q )[(q + M) 2 - m\ 



3E q (q +fi + E q ) 



3E q (q + fi- E q ) 
= 0, 



(A6) 
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where the quantity H(0) (which is nonzero at \x ^ 0) is given in l|27|l. Obviously, the quadratic dispersion law, po ~ p 2 , 
for a massless particle follows from l|A6|l . However, in case of relativistic invariance of the system (i. e. at fi — and 
for coupling constants from the region u l|29|l). or for a color-neutral ground state, when H(0) = 0, the dispersion law 
for a NG-boson changes. Indeed, in this case the term AipoH(Q) from (|A6|) should be replaced by 4ipQH'(0) (H'(0) 
is presented in (H7|)), and one arrives at a linear dispersion law for NG-bosons. 
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